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Abstract. We combine and extend the work of Alexander & Antman 3 and 
Fuller 1111 1121 to give a framework within which precise definitions can be 
given of topological and geometrical quantities characterising the contortion 
of open rods undergoing large deformations under end loading. We use these 
definitions to examine the extension of known results for closed rods to open 
rods. In particular, we formulate the analogue of the celebrated formula Lk = 
Tw + Wr (link equals twist plus writhe) for open rods and propose an end 
rotation, through which the applied end moment does work, in the form of an 
integral over the length of the rod. The results serve to promote the variational 
analysis of boundary-value problems for rods undergoing large deformations. 
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1. Introduction 

In a variational analysis of an elastic structure that is acted upon by end forces 
and moments one needs to consider the work done by the applied loads. The work 
done by the applied force does not usually present any problems. One requires the 
distance travelled by the force, which is usually easy to obtain. More problems 
occur in determining the work done by the applied moment. Here one requires 
the end rotation as 'seen' by the moment, which may lead to complications if large 
deformations are allowed. This paper discusses the ambiguities associated with this 
end rotation and shows how a consistent treatment can be obtained. 

If the structure is very long it may be modelled by an infinitely long rod. This 
has the advantage that powerful techniques from dynamical systems (viewing arc- 
length along the rod as time) [HIE] and variational analysis can be used. The 
natural class of solutions to consider in this case are localised solutions, which de- 
cay sufficiently rapidly towards the ends (other solutions would have infinite strain 
energy and would therefore be non-physical). The boundary conditions to be im- 
posed on such solutions are simple: end tangents are aligned and the ends of the 
rod do not interfere with the localised deformation. The work done by the torsional 
load is then simply the product of applied end twisting moment M and relative 
end rotation R (the angle that has gone into one end of the rod in order to produce 
the deformation, starting from a straight and untwisted reference configuration and 
keeping the other end fixed). 

However, any real-world problem (be it a drill string, a marine pipeline or a 
DNA supercoil) deals with a finite-length rod. For such a rod more complicated 
boundary conditions may be encountered. For instance, the end tangents need not 
be aligned, so that an end rotation that could be used in an energy discussion 
is not straightforward to define. But even in the case of aligned end tangents 
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Figure 1 . Four rod configurations that are not distinguished by a 
prescribed process (homotopy) of boundary conditions from a fixed 
reference configuration: one can go from (a) to (b) either through 
an end rotation or through 'looping'. (After 3 ). 

complications arise if large deformations are allowed. These complications are the 
subject of this paper. 

Let us demonstrate the issue with an example. If we rotate the right end of 
the rod in Fig. OJa) through an angle of — 4tt we obtain the configuration shown 
in (b). If we now move the right end of the rod to the left, the rod pops into a 
looped configuration as shown in (c). If we move the right end further and make the 
loop pass around the right clamp, as in (d), and then pull the ends out again, we 
return to configuration (a) without having rotated the ends. (Note that the process 
illustrated in Fig. UHd) requires that whatever supports the right clamp must be 
released to allow the passage of the rod behind it.) We conclude that we can go 
from configuration (a) to configuration (b) either by end rotation or by 'looping'. 
What, then, is the 'real' end rotation of the deformation (a) — » (b) 'seen' by the 
applied moment and therefore pertinent to an energy analysis? 

It is apparent from this ambiguity that it is not sufficient, a priori, to consider 
only the initial and final configurations of a deformation process to decide on the 
end rotation. We need to be told how the rod was deformed from the one into the 
other, i.e., we need to know the deformation history, not just locally of the ends 
but globally of the entire rod. (This path dependence suggests a relation with a 
geometric phase; details on this are found in |13M181 125] ."I 

Alexander & Antman 3 address the end rotation ambiguity for fixed boundary 
conditions by imagining the open rod to be part of a closed rod. A natural restric- 
tion to a class of open rod deformations is then obtained by demanding that the 
closed continuation should not undergo self-intersections. Applied to Fig. ^ this 
would mean that throwing the rod around the clamp as in (d) is excluded because 
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it requires an intersection with the closure (indicated by the dotted lines). In the 
restricted class of deformations the configurations (b) and (c) can then be assigned 
a unique end rotation of R — —An relative to (a). The ambiguity in the end rotation 
is thus resolved. 

To distinguish between different classes of deformations Alexander & Antman 
use the link of the closed rod, i.e., the topological linking number of two lines 
drawn on opposite sides of the unstressed rod. Since two unknotted closed rods 
with equal linking number may be deformed into each other without undergoing self- 
intersections, there is a one-to-one correspondence between the classes of admissible 
deformations with unique R and the values of the link. The link is thus seen to 
be related to the end rotation. (Note that the —4tt versus ambiguity in the end 
rotation illustrated in Fig.^is reflected in the jump in link by —2 as the closed rod 
intersects itself.) 

The connection between link for closed rods and end rotation for open rods has 
been discussed by many authors. The level of detail varies, but one usually argues 
from the link of a closed rod to the link of an open rod, which is then identified 
with the end rotation. Often this involves the introduction of a closure and the use 
of the celebrated formula ^2 02] 



which expresses the link Lk (of a closed rod) in terms of the twist Tw (a local 
property, in the sense that it can be found by integrating a density along the 
length of the rod) and the writhe Wr. This writhe, which is only a function of the 
centerline of the rod, is not a local property but several expressions exist for the 
calculation of the writhe of an arbitrary closed curve [2] ■ Some of these expressions 
make sense for open curves as well and therefore suggest an extension of writhe 
from closed to open curves |29|. The final result, then, is a formula for the open 
link which in terms of a suitable Euler-angle representation takes the simple form 



This formula cannot be expected to hold true in general since it would make link 
a local property, which it is not. Indeed, the generalisation of writhe from closed 
to open curves used in obtaining is subject to a geometrical condition (see 
Section [2}. Certainly general validity of @ is prevented by the polar singularity 
inevitably associated with Euler angles. In Section 01 formula J2J will be derived 
within a limited class of deformations. 

The present paper improves on previous results in the following ways. 

(1) The closure introduced by Alexander & Antman is effective in the defini- 
tion of admissible deformations if the supports are fixed in both angle and 
position. This is adequate for a large class of experimental situations. If, 
however, the supports are allowed to move, then it is not possible to use one 
and the same closure for all deformations: different deformations require 
different closures. Thus the desired distinction between classes of deforma- 
tions is lost (one can always move the closure along, so that it does not 
'get in the way' during the deformation). We resolve the issue by limiting 
the class of admissible closures in such a way that the separation into dif- 
ferent classes of deformations, characterised by the link, is preserved under 
arbitrary movement of the supports. 

(2) We define a precise class of deformations within which link, twist, and 
writhe arc well-defined. We carefully examine the restrictions imposed on 
the deformations and show in which sense they are necessary. As mentioned 
above, for a consistent definition of link and writhe it is necessary to work 



Lk = Tw + Wr, 
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within a class of homotopies of rods, which connect a given open rod to a 
reference configuration. Despite the requirement of such a connection, the 
newly introduced writhe and link themselves are independent of the choice 
of connection: within the class the writhe and link only depend on the given 
open rod and the reference configuration. 
(3) We show that the link, twist and writhe defined for open rods satisfy the 
classical equality Q. Furthermore, within the class of admissible deforma- 
tions the link is given by (01 ■ In the special case that the end tangents of the 
open rod are aligned the link coincides with the end rotation. Our results 
thus formalise current practice in the literature based on @. However, 
while most applications of formula (J2J can be shown to fall into our class 
of deformations, recent experiments in molecular biology do not always do 
so and care is required in energy discussions (we discuss this in Section 
4). Indeed, we would claim that usage of 10 presupposes the framework 
that we here discuss, and consequently is subject to the limitations that we 
describe. 

The organisation of the paper is as follows. In Section we first define our class 
of open rods. For the elements of this class we define link, twist and writhe and 
derive the extension of Q to open rods. In Section |21 we independently define 
the end rotation and show it to be equal to the open link. We also derive (J2J). In 
Section0]we critically review the defining conditions of the class of rods considered, 
illustrating their relevance with counterexamples. Section [5] discusses our work in 
the light of previous work in the literature. 

2. Results: Link, Twist and Writhe 
For our purposes a rod is a member of the set 

A° = { (r, di) G C 2 ([0, L];R 3 x S 2 ) such that |f| ^ 0, r ■ d x = 0, 

and r is non-self-intersecting } . 

Here and in the following an overdot denotes differentiation with respect to the 
spatial variable s. The curve r is thought of as the centerline of a physical rod (of 
length L) and d\(s) as a material vector in the section at s. As alternatives to 'rod' 
the terms 'ribbon' |121 12"] and 'strip' |3] are also used. A closed rod is an element 
of .A for which begin and end connect smoothly. 

To each point on the centerline of the rod we can attach an orthonormal right- 
handed frame (di(s), d 2 {s), da(s)) of directors by setting 

d a (a) =r(s)/\r(s)\ and d 2 (s) = d a (s) x di(s). (3) 

These directors track the varying orientation of the cross-section of the rod along 
the length of the rod. The twist of a closed rod (r, di) is now defined by 

Tw(r,dx) := — I d 1 (s)-d 2 (s)ds. (4) 
27T J r 

It measures the number of times di revolves around d^ in the direction of d 2 as we 
go around the rod. 

Let r\ and r 2 be two non-intersecting curves. Then the link of r\ and r 2 is 
defined by 

TU \ 1 f j {his) x r 2 (t)} ■ [njs) ~ r 2 (t)} ^ 
The writhe of a closed curve r is 

47r/ r / r \r(s)-r(t)\ 3 
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The argument of this integral is the pullback of the area form on S 2 under the 
Gauss map K 2 — ► S 2 , 



G : (r(s),r(t)) 



r(s) — r(t) 
\r(s)-r(t)\ 



so that the writhe may be interpreted as the signed area on S 2 that is covered 
by this map. For each direction p G S 2 the signed multiplicity of the Gauss map 
(i.e., the number of points (s,i) for which G(s,i) = p, weighted by the sign of 
p- [G s x Gt]) equals the directional writhing number, the number of signed crossings 
of the projection of r onto a plane orthogonal to the vector p ^2G]- m other words, 
the writhe of a closed curve is equal to the directional writhing number averaged 
over all directions of S 2 . 

The link, twist and writhe of a closed rod are related by the well-known Calugareanu- 
White-Fuller Theorem 003111]: 

Theorem 2.1. Let (r,d\) € A be a closed rod as defined above. Then 

Lk(r,d l ) = Tw{r 1 d l ) + Wr(r). (7) 

We review two classical theorems by Fuller which are of interest to us. Note that 
at each point r(s), the unit tangent t(s) = r(s)/\r(s)\ traces out a closed curve on 
S 2 , called the tantrix. Fuller's first theorem relates the writhe of the curve to the 
area A enclosed by the tantrix on S 2 : 

A 

Wr{r) = — - 1 (modd 2). (8) 

Note that the equality modulo two is necessary since the area enclosed by a curve 
on S 2 is only defined modulo Air. 

The second theorem, stated in detail as Theorem 12.61 below, gives under certain 
conditions a formula for the difference in writhe between two closed curves r\ and 
r-i that can be continuously deformed into each another (see Figure |2]): 

Wr( ri )-Wr(r 2 ) = ± J J^- .(i 1+ i 2 ). (9) 



We now proceed with the introduction of the set for which the open link, twist 
and writhe will be defined. For the definition of this set we choose a closed planar 
reference curve r € C 2 ([0, M}; R 3 ) for some M > L. 
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Definition 2.2. 

Al a = { (r,di) G .4° swc/i tftat 3(f, Ji) G C 2 ([0,M] x [0, 1];M 3 x S 2 ) : 

1. /or eac/i A, f(-,A) is an unknotted, non-self-intersecting 

closed curve, 

2. f(s,A) -di(s,A) =0 /orse[0,M], Ae{0,l}, 

3. (f,<Jx)(s,l) = (r,di)(s) /ors6[0,4 

4. f(s,0) = r (s) /or s 6 [0,L], 

5. P(s,0) -r(s,A) > -1 for s G [0,M], A G [0,1], 

6. {f(s, A) : s G [£, M]} is a planar curve for A = and A = 1, 

and these two planes are parallel \. 

Al can be thought of as a class of open rods (r, d\ ) that can be connected by a 
homotopy — satisfying certain requirements — to the reference curve r$. The part of 
the closed rod parametrised by s G [L, M] is called the closure. 

Some of the conditions above are more straightforward than others. Parts 1 and 2 
state that (f , di) is a homotopy of well-behaved closed rods, and by parts 3 and 4 
the homotopy contains the original open rod (r, d\) at A = 1, and the reference 
curve at A — 0. Parts 5 and 6 contain the essential elements of this definition. 
Part 5 is the same non-opposition condition that appears in the statement of Fuller's 
theorem ('Theorem l2.6[) and is required for the conversion of the writhe to a single- 
integral expression. Part 6, which states that the closure should be planar at the 
beginning and the end of the homotopy, is central in the construction. These last 
two conditions are discussed more fully in Section^ 

Note that curves ro exist for which A]. Q is empty: if the three vectors ro(0), 
fo(L), and ro(L) — ro(0) are independent, then the open curve ro can not be closed 
by a planar closure, so that the set of homotopies with planar closures that connect 
to ro is empty. 

We are now in a position to define the new functionals open link, open twist and 
open writhe for open rods. 

Definition 2.3. Let (r, d\) be a rod in A\ - Then the open twist of (r,di) is 

L 



Tw°{r,dx) := / di ■ (f x d x ) ds, 
27T Jo 



the open writhe of r is 
Wr°(r) := Wr (f( , 1), = J- f" f" [nM)-nU)]-^.l) xf(M)] ^ 



i*J« J„ |f(»,l)-f(t,l)P 
and the open link of (r, di) is 



Lk°[r,d{) :=Lk{f{;l),d 1 {-,l))-^- / d x {s, 1) • (f (s, 1) x dt(s, 1)) ds. 

Zn 



M 

L 



Note that in the last definition we subtract any twist the closure might have. It 
follows directly from the construction that the new concepts also satisfy the classical 
relationship: 
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Corollary 2.4. Let (r,di) G A]. Q be an open rod. Then 

Lk°(r, di) = Tw°{r, di) + Wr°(r). 

Theorem 2.5. For any open rod (r,di) G A], , the open twist, writhe, and link are 
well-defined. 

For writhe and link this is non-trivial, as different homotopy closures (f, d\) 
might be expected to give rise to different values. 

Proof. We first state Fuller's second theorem in a more precise form. 

Theorem 2.6. Let r\ (0 < A < 1) be a homotopy of closed non-self-intersecting 
curves, regularly parametrized with a common parameter s G [0, L] . Let t\ be the 
tantrix of r\. Ifto(s) • t\(s) > —1 for all s G [0, L], A G [0, 1], then 

Wr(n) - Wr(r ) = ±- t *° * h ■ (i + h)ds. (10) 
Z7r Jo 1 + £o ■ £i 

To our knowledge, a rigorous proof was first given by Aldinger et al. [2]. 

To prove Theorem 12.51 for the open writhe, let (f, rfi) be a homotopy associated 
to (r, di). By definition, f(-,0) and f(-, 1) are planar for s G [L,M]. Let us denote 
the planes by Vq and Vi\ these are parallel by Definition 12. 21 6. Let V be the plane 
through the origin parallel to both. 

We have defined the class of open rods A\ such that Theorem l2.6l can be applied. 
Denote the tantrices of f (•, 0) and f(-, 1) by to and t\ respectively. Then 

Wr(f(;l))-Wr(f(;0)) = 1 f M to(s) x Ms) + . & 

27r J 1 + t (s) • ii(s) 

The argument of the integral vanishes for s G [Z, M]: since to(s), fi(s) G 1^ for 
s e [L, M] we have t (s) x ^(s) 1 V and t (s) + ii(s) G V for s G [L, M}. Hence 

[*o(«) x*i(*)] ' [i (s)+ii(s)] =0. 
Moreover, since f(-,0) is planar, H / r(f(-,0)) = 0. We conclude 

Wr(f(.,l)) = * f L /° (a) x tl(a) ■(i ( a ) + i 1 ( a ))d a . (12) 

Since this integral only depends on the reference curve and the open rod itself, and 
is otherwise independent of the choice of closure and homotopy, this proves the 
claim for the writhe. 

For the link, let (f, d\) be a homotopy associated to (r, d\) by Definition 12.21 
Denote the closed curves f(-, 1) and Ji(-, 1) by f and d\ respectively. We denote 
the open twist evaluated over an interval s G [a, b] by Tw° a Then 

Tw °L,M](r,di) = Tw °oM]( r 'di) -Tw\ QiL] (r,d x ) 
= Tw [o,M]( r ,di) -Tw°{r,d{). 

Hence 

L*°(r,di) = LKMi)-Twf L!M] (M!) 

= Lk{f,d 1 )-Tw{f,d 1 )+Tw°{r,di) 
= Wr(f) +Tw°(r,d 1 ) 

= Wr°{r) +Tw°(r,d 1 ). (13) 
It follows that Lk°(r,di) is independent of the chosen closure (f,d\). □ 
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Equation (I12|) is an important motivation of this work, since it expresses the 
writhe in terms of a single rather than a double integral. For the purpose of 
variational analysis this is an obvious advantage. It is especially useful when the 
link, and therefore indirectly the writhe, can be identified with the rotation of the 
ends; this requires that the end tangents remain equal throughout the deformation, 
and this case is treated in the next section. 

Remark 2.7. For simplicity we have chosen to introduce one class A 1 as the basis 
for the definitions of open link, twist and writhe. For each of the three definitions 
separately, however, not all of Definition 12.21 is required. Open twist can be defined 
directly in terms of (r, d\), without the need for an extension; open link requires a 
closure, but no homotopy; and open writhe requires a homotopy, but the director d\ 
can be disposed of. 

Remark 2.8. A natural question to ask is whether the open link, twist and writhe 
reduce to their classical counterparts when an open rod is transformed into a closed 
rod by lining up and connecting the ends. This is not the case, as we demonstrate 
in Section 0] 

3. Results: End-rotation and Euler angles 

It is common in applications to assume that the end tangents of the buckled rod 
are kept constant and equal during the deformation process. For comparison with 
an end rotation we introduce this additional condition. Throughout this section 
we also assume that ro|[o,L] is straight; without loss of generality we assume that 
ro is a constant unit vector v on [0, L]. Finally, again without loss of generality we 
choose the director d\ constant on the reference curve ro|[o,L] : 

Definition 3.1. 

A 2 ro = { Mr) £ Al : ro(<0 = v e S 2 for all s £ [0,L], 

r(0, A) = f(L, A) = r(M, A) = v for all A £ [0, 1], 
Ji(s,0) = Ji(0,0) for all s £ [0,L] }. 
The following formula is a direct consequence of Theorem 12. 61 
Corollary 3.2. Let (r,d\) £ A 2 , and let t be the tantrix of r. Then 

Wr°(r) = ± [ \ VX ^l .i(s)ds. (14) 

27T J Q 1+Vt(s) 

In the present case of a straight rolron the dependence on ro of the open writhe 
of a given open rod takes a particularly simple form: 

Theorem 3.3. Under the conditions of Corollary \3. 21 let £1 = S 2 \ {— t(s) : s £ 
[0, L]}. Then the function 

I [ L vx t(s) ., , , 

veS 2 * — > — / V- • t(s) ds 

2irJ 1 + v t(s) V ' 

is constant on connected components of ft. 

The proof is given in the appendix. The interpretation of this theorem is as 
follows: when the end tangents are aligned, the tantrix given by the rod (without 
closure) forms a closed curve on S 2 . The integral above represents 'area enclosed 
by the curve' for a given 'choice of area' (cf. (jHJ). When the vector v crosses the set 
{— t(s) : s £ [0, L}} the geodesic connections between v and t(s) change direction, 
causing the integral to represent a different choice of area, and therefore causing 
the integral to jump by 4ir. 



LINK AND WRITHE FOR OPEN RODS 



9 



With fixed end tangents we can introduce a fourth quantity, the end rotation. 
We denote d(-)/d\ by d\(-). 

Definition 3.4. Let (r,di) G A?j. , and let d,3(-,-) — r(>, •)/ -)|, J2 = 4 x d\- 
We define the end rotation by 

R(r,di):= f d\di{L,X)-d 2 {L,\)d\- [ d\dx(0, A) • d 2 (0, A) d\. 
Jo Jo 

To study the relationship between end rotation and open link, twist and writhe 
we introduce a particular choice of Euler angles for an open rod (r, d\ ) . Recall that 
for every s G [0, L] there is an orthonormal director frame (di(s), d 2 (s), °fa( s ))- We 
express this frame in terms of angles 9, ip, 4> with respect to a fixed basis (ei, e 2 , 63) 
as follows 

d\ = {— sin tp sin <p + cos tp cos <f) cos 6) e%+ 

(cos tp sin <f) + sin -0 cos (p cos 9) e 2 — cos <p sin 9 e3, 

d 2 = (— sin tp sin — cos sin </> cos 9) e\+ (15) 
(cos ip cos + sin ip sin cos 9) e 2 — sin </> sin e3 , 

g?3 = cos ^ sin 9 e\ + sin ^ sin e2 + cos 9 . 



This choice of Euler angles follows Love JT2, art. 253]. For rods in the class A^. 
we choose e^ parallel to v\ note that by this choice the non-opposition condition 5 
in Definition 12.21 coincides with avoidance of the Euler-angle singularity at 9 = n. 
Therefore the smoothness assumptions on (f, d\) in A^ imply C^-regularity for (p, 
ip, and 9. 

Lemma 3.5. Let (r, d\) G A 2 TQ he an open rod with an associated homotopy (f, d\), 
and let d 2 and d% be constructed from f and d\ according to |3J). Let (p,ip,9 : 
[0, M] x [0,1] — > M be the Euler-angle representation of {di 1 d 2 ,d^). Then 



R{r,dx) = I [4>(s,l)+ip(s,l)]ds. 



From this Lemma we conclude 



Corollary 3.6. For open rods (r,dx) G A^ , R(r,d%) is independent of the choice 
of extension (f,di). 

Proof of Lemma \S.5\ Using the definitions of the Euler angles, we find 

d\d~i(s, A) • d 2 (s, A) = d\4>(s, A) + d\ip(s, A) cos#(s, A). 

For s = 0, L and A G [0,1] we have set d 3 (s, A) = e 3 , and hence 9(s, A) = for 
s = 0, L; therefore d\d~i(s, A) ■ d 2 (s, A) = d\<p(s, A) + d\ip(s, A) for s = Q,L. 

Since + %p is a continuously differentiable function on V := [0, L] x [0, 1], the 
integral of the tangential derivative of cp + ip along dV vanishes: 

/ |-(0 + VO = o, 

JdV OT 

where r is the clockwise-pointing unit vector tangential to dV . 
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Hence 

[d s <t>{s, 1) + d a ip( S) 1)] ds- [d s cj>(s, 0) + d a ip(s, 0)] ds = 
n Jo 

[d x <j>(L,\) + dxi>(L,A)]d\- f [d x <l>(0,\) + dx4>(0,X)]d\ 
i) Jo 
i S=L 

[a A 0( s ,A) + 9 A V(s,A)]dA 



o 



s=0 



d\di(s, A) • d 2 (s, A) dX = R(r, d%) 
=o 



□ 



Within this framework Euler-angle formulae for open link, twist and writhe are 
obtained: 

Lemma 3.7. Let (r,di) be an open rod in A% Q . Then 

Wr°(r) = ^- [ ip(s)(l-cos6(s))ds, (16) 

and 



Tw°(r,di) = ^- [ [<j)(s)+ip(s)cos6(s)}ds. (17) 
2tt Jo 

Proof. The formula for twist is easily found by using (|15|l in the definition of twist, 
as in the proof of Lemma 13.51 For the writhe we apply Corollary 13.21 and use the 
fact that v — e$. □ 

The main result of this section states that for rods in A% end rotation is equal 
to the open link: 

Theorem 3.8. Let (r,d±) be an open rod in A^. Q . Then 

R{r 7 di) = 2nLk {r, d x ) = f [0(s) + rp(s)] ds = 0(L) + ijj{L) - 0(0) - ^(0). (18) 

Jo 

Proof. By Corollary 12 . 41 and Lemma 13. 71 we obtain 

Lk°(r, di) = Wr°(r) + Tw°(r, d x ) = ^- [ [<j)(s) + ip(s)] ds. 

27r Jo 

Since Jq L [0(s) + ip(s)] ds = R by Lemma f3. 51 we have the desired result. □ 

4. Critique of the approach 

The example of Fig. shows that end rotation can only be defined for a defor- 
mation history. For the purpose of analysis of elastic structures this dependence on 
deformation history is undesirable. The approach in this paper, which is shared by 
many others (see the next section), is therefore to construct a class of deformation 
histories (homotopies) within which the end rotation can be expressed in terms of 
the initial and final states only. In this section we critically review the essential 
ingredients of this approach. 

The closure. We obtain a separation into deformation classes by the introduc- 
tion of a closure. The classes are characterised by the link of the closed rod-closure 
combination. The price we pay with this construction is the dependence on the 
choice of closure, which at first sight might seem to be a defect of the formula- 
tion. As Alexander & Antman |3] point out, however, this dependence is entirely 



LINK AND WRITHE FOR OPEN RODS 



11 



natural: the precise form of the closure can be regarded as describing the way the 
rod is supported. Different systems of supports necessarily allow different classes 
of deformations. 

Although in most applications throwing the rod around the clamp as illustrated 
in Fig. \TJ[d) is physically prevented, there do exist exceptions to this rule. In some 
recent experiments DNA molecules are manipulated with the help of a magnetic 
bead attached to the end of the molecule and held in a magnetic trap which allows 
the simultaneous application of a force and a moment |3()| . In this case no rigid 
mechanical support is present and the molecule is free to loop around the beaded 
end. If this happens repeatedly, then the end can rotate over an arbitrarily large 
angle under the applied moment without a concomitant change in configuration 
(cf. Fig. ^ after a rotation of the end by —4tt the rod returns to its original 
shape). Thus one might estimate the wrong energy, by ±47rM for each 'looping', if 
a deformation would go outside the class, i.e., if link was not conserved. However, 
Rossetto & Maggs in a recent paper show that for micron-sized beads the 
applied tension in many experiments is large enough (on the order of femtonewtons) 
to make these link violations rare, and they proceed to introduce a closure to study 
configurations of constant link. 

Incidentally, the fact that a rotation of ±47r, and not one of ±27r, brings one back 
to where one started has its origin in the topological nature of SO (3), the group 
of rotations in R 3 . Specifically, SO(3) is not simply connected: for every rotation 
R G SO (3) there are two homotopy classes of paths from the identity of the group 
to R. This means that for a given rod orientation there are two distinct classes of 
configurations for the rod which cannot be deformed into each other while keeping 
the ends fixed. Rods with any even number of end turns (including zero) lie in one 
class; rods with any odd number of turns lie in the other. The same topological 
property forms the basis of the famous Dirac Belt Trick JSJ, which in classrooms 
is often illustrated by rotating a cup, held in the palm of one's hand, twice around 
a vertical axis by a suitable motion of the arm to bring both cup (with contents) 
and arm back to their initial positions (see ^01 for a demonstration). 

The reference curve. The class A], is defined for a fixed closed reference curve 
ro (which may or may not include the unstressed centerline of the rod). The open 
writhe Wr° will, in general, depend on the choice of this curve; on the one hand, 
by the fact that the reference curve restricts the class of admissible homotopies via 
the non-opposition condition, and on the other hand, by the explicit dependence 
on t in l|12|l. Similarly, the end rotation R will depend on r , as is to be expected 
since R is defined (in Definition 13. 4|) as the end rotation incurred in deforming 
(r ,<ii(-,Q))|[o,L] into (r,di). 

For certain cases, however, the dependence can be described more precisely, as in 
Theorem 13. 31 for straight reference curves. Though beyond the scope of this paper, 
it is also possible to prove a similar result under rotation of more general reference 
curves. For instance, one could extend Definition ^. 41 Lemma and Corollarv l3.6l 
to a larger class than A% by requiring of ro|[o,L] only that its end tangents be 
equal. A complication would arise, however, in that violation of the non-opposition 
condition would no longer coincide with the Euler-angle singularity. Consequently, 
the non-opposition condition would no longer assure us of ^-regularity of cj>, ip, 
and 9. 

The non-opposition condition cannot be dropped. The non-opposition 
condition listed in the definition of A\, Q (condition 5) is imposed by the application 
of Theorem l2.6l This might seem to be merely a technical restriction: after all, the 
open writhe is defined in terms of the writhe of the closed curve fDefinition l2.3[l and 
the latter is well-defined even if, somewhere along the homotopy, the non-opposition 
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condition is violated. Therefore it might be expected that the statement of well- 
posedness holds true without condition 6 (even though our proof evidently does 
not), and that only non-self-intersection of the closed structure is required. 

In fact the situation is not that simple. Figure shows homotopy paths con- 
necting the reference configuration (a) with the deformed rod-closure combinations 
(b), (c) and (d), where the rod itself (represented by the thick line) is the same in 
each of the three deformed states. We can imagine the deformed rod to be nearly 
planar, with the two strands crossing at a short distance from each other. Then 
the rod and its closure in case (b) have writhe close to — 1. In case (c) one adds 
or subtracts 1 to the writhe of the rod-closure combination for each full turn of the 
end. The writhe of the combination can therefore be made arbitrarily large. In 
case (d), finally, the writhe is close to —2. 




Figure 3. An example showing that the non-opposition condi- 
tion 5 in Definition 12.21 cannot be disposed of. The three final 
states (b), (c) and (d), with identical shapes for the open rod, 
have different values of writhe. In going from (b) to either (c) or 
(d) the non-opposition condition is violated. 



It is not difficult to see that one may construct a homotopy between (a) and (b) 
without violating the non-opposition condition, provided the loop has been twisted 
through an angle strictly less than tt. Since the continuation homotopies to (c) and 
(d) satisfy all conditions of Definition 12 . 21 other than the non-opposition condition, 
it follows from Theorem 12 . 51 that these homotopies cannot be constructed without 
violating the non-opposition condition. This may also be verified by inspection. 

This example shows that simply removing condition 5 from Definition 12 . 21 leads 
to ambiguities in the definition of writhe (and therefore of link). The example also 
suggests that if a well-defined writhe is to be constructed without the inclusion of 
the non-opposition condition, then additional restrictions must be imposed on the 
closure. In homotopy (c) the closure remains planar throughout the homotopy, but 
the end tangents vary; in homotopy (d) the end tangents are constant, but the 
closure is only planar at the beginning and the end of the homotopy. To rule out 
homotopies (c) and (d) (necessary for a well-defined writhe) we can require the end 



This may be verified by using the characterisation of writhe as the average of the directional 
writhing number, as explained in Section |2] This number is determined by counting signed 
crossings in a projection of the curve onto a plane. 
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tangents to be fixed and the closure to be planar throughout the homotopy. It is 
possible that for a well-defined writhe further conditions must be imposed. 

Euler-angle singularity vs. the non-opposition condition. When the ref- 
erence configuration is straight and end tangents remain constant during the ho- 
motopy, the non-opposition condition is equivalent to avoidance of the Euler-angle 
singularity at 9 — it. Although this is partly a coincidence, the two issues both 
stem from the topological properties of S 2 . 

The Euler-angle singularity results from the fact that S 2 is not homeomorphic 
to (any part of) R 2 . Any parametrization of S 2 by a single cartesian coordinate 
system will therefore have at least one singular point. On the other hand, the 
non-opposition condition is necessary — in this article — for the single-integral rep- 
resentation of writhe of Theorem 12. 61 In this representation the ambiguity of area 
'enclosed' by a curve on S 2 is resolved by taking a perturbation approach. The non- 
opposition condition is the realization of the unavoidable limits of this approach, 
and therefore again stems from the topology of S 2 . 

As mentioned above, however, the non-opposition condition remains an unsatis- 
factory element in the definition of open writhe. Perhaps a concept of open writhe 
is possible that bypasses this condition. 

Open writhe is not rotation invariant. The definition of A\ depends on 
the choice of the reference configuration. For a given reference configuration, an 
open rod in A\ may not be freely rotated without leaving A]. Q . This is readily 
demonstrated by rotating the reference configuration itself: after a rotation of it 
about an axis perpendicular to the plane of the reference curve the non-opposition 
condition is violated at every point on the curve. 

This may lead to surprising results. In Figure0]two homotopies are shown. The 
first is a variation on homotopy (b) of Figure |3 while in the second we lengthen 
the open-rod part and shorten the closure part. In addition, we construct the 
homotopies such that the final configurations are close, up to a rotation (emphasized 
by the mark at one end of the open rod). In (a) the open writhe is close to 1, while 
in (b) it is close to 0. 



Figure 4. Two elements of A), that differ only by a rotation, but 
for which the writhe is different. The dot emphasizes the difference 
in orientation. 

This remark also resolves an issue raised in Remark 12.81 does the open writhe 
change continuously into the classical writhe for closed rods, when an open rod is 
transformed into a closed rod by lining up and connecting the ends? The answer is 
no — for the resulting closed writhe would be rotation- invariant, contradicting the 
remark above. 
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5. Discussion 

The topological issues associated with large deformation discussed in this paper 
are not of great concern in more traditional engineering applications. As long as de- 
formations are such that the integral in @ remains well-defined for a suitable choice 
of Eulcr angles (i.e., as long as the angles stay away from the polar singularity) open 
link and end rotation are given by (0). However, in more modern applications of 
structural mechanics, such as in molecular biology, large deformations occur more 
routinely and more care is required. Indeed, Fuller's 1971 paper jll| was inspired 
by supercoiling DNA molecules. In this paper the author also already introduces a 
(planar) closing curve in order to compute the writhe of a simple (infinitely long) 
helix. Following the pioneering work of Fuller an extensive literature has emerged 
on the application of elastic rod theory to DNA supercoiling (c/. the survey article 
by Schlick gHj). 

Open rods have become popular models for DNA molecules since by the early 
1990s single-molecule experiments have become possible. First this involved an 
applied force only [J7J; later, once the molecule could be prevented from swivelling 
at its (magnetically) loaded end, this involved both an applied force and an applied 
moment j2D]. Analytical studies have addressed DNA in isolation as well as 

in dilute solution using a statistical mechanics approach 0H1 1201 H] • 

Bcnham 4 appears to have been be the first one to write down an isoperimetric 
variational problem based on J3J) in order to find equilibrium configurations subject 
to constant link. He first uses Fuller's result JJjJ with r-i taken to be a suitable closed 
planar curve. Then Wr{r%) = 0, and one obtains a single-integral formula for the 
writhe of the curve in question r% . When combined with the single integral for the 
twist (@J), this leads, via 0J, to a single- integral expression for the link, which in a 
suitable Euler- angle representation is given by J2J- Bcnham then observes that "the 
integral expressions for Lk, Tw and Wr may be constructed regardless of whether 
the structure is closed". Many workers have since followed Benham's example to 
write link and writhe of an open rod as single integrals. 

The implicit assumption in this approach is that there exists a continuous de- 
formation from r\ to the planar curve which avoids opposition of corresponding 
tangents. This may not be obvious, since the non-opposition condition must be 
applied to the closed combination of rod and suitable closure. In Section0]we have 
given examples of what can go wrong if the condition is not satisfied. In addition, 
the example of Figure 0] shows that even within the limits of the non-opposition 
condition the writhe is not invariant under rotation, implying that this approach 
may lead to counterintuitive results. 

Frequently, the end rotation needed in an energy analysis is simply assumed to 
be equal to the link as given by J2J (e.g., [HI GO!)- We have shown that end rotation 
can only be defined in a consistent way within a class of homotopies of rods. Such a 
class is constructed with the help of a closure. We define end rotation independent 
of link and show the two to be equal within a suitable class of allowed deformations. 
We should also remark that our closure is a rod (r, d±) rather than just a centerline r. 
Most authors initially only introduce a closed centerline, which makes the writhe 
well-defined, and subsequently assume the closure to be twistless if a link or end 
rotation is required. We formalise this practice by explicitly specifying a d\ for the 
entire closed structure. 

Knowing the precise restrictions on allowed deformations is important in statisti- 
cal mechanics studies. To obtain the correct averages one must consider ensembles 
of admissible configurations. In numerical computations this means that one must 
take care to simulate configurations (through a Monte Carlo algorithm, a 'growth' 
algorithm or otherwise) with the right topological constraints. Specifically, one 
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wants the configurations to be unknotted and to have constant link (although it is 
good to remember that DNA in its natural environment functions in the presence 
of topology-changing enzymes). This means that one must forbid self-crossings of 
the configuration as well as crossings of an (imaginary) closure. Mindful of this, 
the authors in graft the ends of the molecule to an external surface and run 
'sticks' out from the ends of the molecule to infinity, thereby 'virtually closing' the 
generated chain at infinity. A similar construction is used in |25| . In this latter 
work knotted configurations are not eliminated, it being argued that the molecular 
statistics is dominated by unknotted configurations. 

In order to avoid the awkward non-opposition condition in Fuller's second the- 
orem there have been direct approaches via the double integral © instead. In 
|291 HT\ simple shapes are considered with planar closures for which the integral 
can be evaluated explicitly. It is then shown that the closure gives a relative contri- 
bution to the writhe which tends to zero as the length of the rod tends to infinity. 
The double integral is also used in the numerical study in where it is shown 
that the contribution to Wr from the interaction of the closure with the basic chain 
is of the order of 1%. Various numerical schemes for the computation of the writhe 
double integral for a discretised curve are discussed and compared in |lf>| . Useful 
rigorous error bounds on numerically computed values of Wr based on polygonal 
(i.e., piece- wise linear) approximation are given in [7|. 

Our approach to a consistent definition of link, writhe and end rotation is firmly 
based on the introduction of a closure. There have been various formulations of 
writhe for open curves without the use of a closure. One approach, especially taken 
in knot theory and in studies of self- avoiding chains, is based on the characterisation 
of writhe as the average over all planar projections of the sum of signed crossings JT] 
(e.g., 12210)- This approach does not require a closure and can be applied to curves 
with arbitrary end tangents. It could form the basis of an alternative extension 
of Q to open rods. (Here we can remark that link also has an interpretation in 
terms of crossing numbers, namely: the linking number of two curves is equal to 
the number of all signed crossings of the curves in a regular planar projection, i.e., 
one satisfying a transversality condition; see [211 ■) The precise relation between the 
writhe of an open curve obtained via these planar projections and the open writhe 
obtained by using the tantrix area on S 2 , or the double integral ©, is still an open 
problem. 

If an exact writhe is not required and the fractional part modulo 2 is sufficient 
then formula ijHJ) in terms of the area enclosed by the tantrix on the unit sphere 
can be used. Rossetto & Maggs [25] point out that this approach can also be 
used to generalise the writhe to curves whose end tangents are not aligned and 
therefore have open tantrices. By exploiting the connection between writhe and a 
geometric phase they show that the canonical way to close the curve is by means 
of a geodesic (great circle). This geodesic is unique as long as the two end points 
are not antipodal. The fractional writhe is then again given by the enclosed area. 
This prescription is used by Starostin to derive results for the writhe of smooth 
as well as polygonal curves. Cantarella |7] generalises Fuller's second theorem, and 
with it the spherical area formula (JHJ), to polygonal curves. 

Appendix A. Proof of Theorem 13.31 

Let / : S 2 — > K be the function mentioned in the assertion. Pick vo £ ft and let 
51q be the connected component of containing vq. Define the set 



A = {v £ n n : f(v) = f(v )}. 
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The function / is continuous on implying that the set A is relatively closed in 
Slo- We will show below that / is constant on all open balls B C VLq, implying that 
A is also open. Since A is non-empty it follows that A = ilo and the Lemma is 
proved. 

For a given vector w £ S 2 , let denote the rotation around to through an 
angle <fi. We fix the direction of rotation in the following way: with respect to an 
orthonormal basis (u>, w,lo x w) for a suitable w € S 2 , write as 

/I \ 

R4, = cos cj> — sin . 
V sin0 cos0 J 

With this choice, 

— i?0f|0 = o = oj x v for any v £ S 2 . (19) 

Set v<j> = R<pv. Using equation (|19() . we have 

d v^xt ■_ [(1+^- t)((u) x v$) x t) - x t)((ui x ^) ■ t)} . 
dej) 1 + Vrj > ■ t ' (1 + v^-t) 2 ' ' [ ' 

Setting 7 = \lo x we introduce an orthonormal coordinate system 

ei=v ( f > , e 2 =7 _1 wx^, e 3 = 7" 1 ^x(wxu ( i), 

and we write tx, ti, H for the coordinates of t with respect to this basis; these are 
functions of the curve parameter s. The right-hand side of 1)20(1 becomes 

(1 + tiXjaji - jijg) - t 2 (t 2 4 - ^^2) 

7 (l + ii) 2 

Using the equalities t 2 + t\ + t 2 = 1 and tit\ + t-iti + ^3 = this is seen to be 
equal to 

d t 3 
7 ds l+tx' 

Therefore 

d<p J l + v^-t J ds l + tx 

The last equality results from the assumption of aligned end tangents. This proves 
the Theorem. 
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